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ABSTRACT 

Since the launch of the Fermi 7 - ray Space Telescope on June 11, 2008, the LAT instrument has 
solidly detected more than 20 GRBs with high energy photon emission above 100 MeV. Using the 
matched filter technique, 3 more GRBs have also shown evidence of correlation with high energy pho- 
ton emission as demonstrated by Akerlof et al. In this paper, we present another GRB unambiguously 
detected by the matched filter technique, GRB 081006A. This event is associated with more than 13 
high energy photons above 100 MeV. The likelihood analysis code provided by the Fermi Science 
Support Center (FSSC) generated an independent verification of this detection by comparison of the 
Test Statistics (TS) value with similar calculations for random LAT data fields. We have performed 
detailed temporal and spectral analysis of photons from 8 keV up to 0.8 GeV from the GBM and 
the LAT. The properties of GRB 081006A can be compared to the other two long duration GRBs 
detected at similar significance, GRB 080825C and GRB 090217A. We find that GRB 081006A is 
more similar to GRB 080825C with comparable appearances of late high energy photon emission. 
As demonstrated previously, there appears to be a surprising dearth of faint LAT GRBs, with only 
one additional GRB identified in a sample of 74 GRBs. In this unique period when both Swift and 
Fermi are operational, there is some urgency to explore this aspect of GRBs as fully as possible. 

Subject headings: gamma rays: bursts 



1. INTRODUCTION 

Gamma-Ray Bursts (GRBs) are extremely luminous 
explosions in the Universe but most observations are lim- 
ited to photons below a few MeV. High energy emission 
above 100 MeV was only detected a few times by the 
EGRET instrument (Dingus ct al. 1995) and now more 
recently by AGILE (Giuliani et al. 2008). The recently 
launched Fermi 7 - ray Space Telescope has enlarged the 
opportunity to study high energy radiation from GRBs 
by providing large apertures and wide fields of view. 
Two onboard instruments, the Gamma-ray Burst Mon- 
itor (GBM; Meegan et al. 2009) and the Large Area 
Telescope (LAT; Atwood et al. 2009), overlap energy 
bands to span from 8 KeV to above 100 GcV. With other 
satellites providing more precise localization ability (e.g. 
Swift; Gehrels et al. 2004, INTEGRAL; Winkler et al. 
2003, AGILE; Giuliani et al. 2008), we have a unique 
opportunity to study the physical mechanisms of GRBs 
across a very wide dynamic energy range. 

The Fermi /LAT covers the energy range from below 
20 MeV to more than 300 GeV with an effective field-of- 
view (FOV) - 2.4 sr (Atwood et al. 2009). The other 
Fermi instrument, the GBM, is sensitive to the 8 keV 
- 40 MeV range and covers the entire unocculted sky. 
The GBM identifies GRBs in real time (Meegan et al. 
2009) with a rate of ~ 250 events per year (Paciesas et 
al. 2010). During its first 27 months operation, ~150 
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GRBs have been detected by the GBM at locations that 
were simultaneously within the Fermi /LAT field of view. 
However, of these simultaneously observed events 5 , only 
20 have been detected at a threshold of more than ~ 
10 high energy photons above 100 MeV (Granot 2010), 
corresponding to a rate of ~ 9 GRBs per year. 

Using the matched filter technique, the Fermi de- 
tection threshold was reduced to a level of ~ 6 high 
energy photons with the concomitant identification of 
3 additional detections as reported by Akerlof et al. 
(2010, 2011, hereafter A10, All), namely GRB 080905A, 
091208B and GRB 090228A. In Section 2, we present the 
first high energy photon detection of GRB 08 1006 A using 
this matched filter technique. The detailed observations 
obtained by the GBM and LAT, including temporal and 
spectral properties, are described in Section 3. The com- 
parison of GRB 08 1006 A with other two GRBs detected 
with similar fluences are discussed in Section 4 and sum- 
marized in Section 5. 

2. DETECTION OF GRB 081006A USING THE MATCHED 
FILTER TECHNIQUE 

Following the method developed by Akerlof et al. in 
A10 and All, we apply the matched filter technique to 
an extended GBM sample by lowering the GBM fluence 
threshold below the 5 /ierg/cm 2 minimum that defined 
the GRBs listed in Table 1 of A10. The extended GBM 
sample consists of cither GRBs that have GBM fluences 
less than 5 //erg/cm 2 or have no GBM fluence informa- 
tion whatsoever. The remaining criteria are the same 
as in A10. This new sample of 74 GRBs, listed in Ta- 
ble 1, was searched for high energy photon emission us- 
ing the matched filter technique. One GRB, namely 

5 http: / /fcrmi. gsfc.nasa.gov / ssc / observations / types / grbs /grb_tablc/ 
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GRB 081006A, was revealed with a large matched weight 
value, (J2 w i = 652.460, more than 26 times larger than 
any other GRBs in this set. The statistical significance 
of this detection can be estimated from the cumulative 
distribution of matched filter weights for a large sample 
of otherwise similar LAT fields. As shown in Figure 4, 
after including the GBM-triggcred sample size of 74, the 
probability of such an event arising from random back- 
ground photons is significantly less than 0.4%. A more 
detailed discussion is provided in ^3.21 

The triplet cluster finder associated with the matched 
filter technique estimates the location of GRB 081006A 
to be RA = 134.4° and Dec = -61.8° (J2000.0) with an 
uncertainty of ~ 0.5°. 13 high energy photons are spa- 
tially clustered in the nominal GRB 081006A direction. 
As shown in Table 2, 11 of these photons lie within a 
2° radius. Figure 1 shows the LAT photon sky map of 
GRB 081006A. Our estimated direction (filled blue dot) 
is 6.6° away from the GBM determination (center of the 
map) which has ale uncertainty of 4.5° with an ad- 
ditional systematic uncertainty of the order of 2.5° (van 
der Horst 2008), thus demonstrating consistency at the 
1 a level. We note that GRB 081006A occurred -120s 
before the LAT ceased to take data as the satellite moved 
closer to the South Atlantic Anomaly (SAA). This must 
be accompanied by an increased photon background rate 
but over the whole LAT field before and after the GRB 
trigger we find no significant rate change that would cre- 
ate a spurious event that could mimic the signature of 
GRB 081006A. 




Fig. 1.— LAT high energy photons sky map for GRB 081006A. 
The diameter of each dot is proportional to its statistical weight. 
Thus, the largest diameters represent Event Class 3, etc. The 
dotted circles around each point indicate the 1— ct errors. The figure 
is centered on the nominal coordinates furnished by the GBM; the 
blue dot on the lower left shows the GRB coordinates computed 
by the cluster algorithm described in the text. The large green 
circle depicts the boundaries of the 16.0° radius cone that defines 
the fiducial boundaries for the cluster search. The plot axes are 
aligned so that North is up and East is to the right. 
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List of 74 GBM triggered GRBs 
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TABLE 2 

List of high energy photons for GRB 081006A 
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a photon ID number 
b time after the trigger 

c distance to the new estimated GRB location 

d energy of each photon 

e photon class 

f weight of each photon 

indicates diminished we for highest energy triplet cluster pho- 
ton 

3. FERMI GBM AND LAT OBSERVATIONS 

3.1. GBM observations 

At 14:29:34 UT on 6 October 2008, the Fermi 
GBM triggered and located GRB 081006A (trigger 
244996175/081006604; van der Horst 2008). Using the 
GBM trigger data, the on-ground calculated location was 
RA = 142.4°, Dec = -67.4° (J2000) with an uncertainty 
of 4.5° (1 a containment, statistical only). The GRB Tgo 
is ^7 s long (van der Horst 2008), clearly detected by Nal 
detector numbers and 3. The BGO detector located on 
the same side, B0, also detected the increased flux. 

Figure 2 shows the light curves obtained by the GBM 
and LAT. The upper plot is the background-subtracted 
light curve for the two Nal detectors, NO + N3, in the 8 
keV to 800 keV energy range. The counts were binned 
in 0.5 s intervals to provide a good visual representa- 
tion of the intensity behavior. The middle plot is the 
background-subtracted light curve for the BGO (B0) de- 
tector in 280 keV to 4 MeV energy range with the same 
time binning. The GBM data shows a single peak con- 
centrated in the first 2 s after the trigger. The lower plot 
is the light curve for the LAT detector above 100 MeV 
within an angle of 12° with respect to the nominal GBM 
direction described above. 

3.2. LAT Observations 

The boresight angle of GRB 081006A with respect to 
the LAT was 16° at the time of the trigger. The LAT did 
not slew to the burst direction but remained in normal 
observation mode, taking data for ~120 s before entering 
the SAA. GRB 081006A was not reported by the LAT 
team but we have shown that it is unambiguously iden- 
tified by the matched filter technique. With at least 13 
high energy photons detected, it is bright enought to al- 
low a likelihood analysis using the standard Science Tools 
software package 6 provided by the FSSC. 



The gtlike tool, part of the standard Science Tools 
software package, generates the Test Statistic TS = 
2Alog(likclihood) to compare models with and without 
an assumed source. The TS value associated with each 
source is a measure of the source significance or cquiv- 
alcntly the probability that such an excess can be ob- 
tained from background fluctuations alone. A TS value 
of 25 corresponds to a 4.6cr significance, approximately 
the square root of TS (Abdo et al. 2010a). 

The first step to perform a likelihood analysis with the 
gtlike tool is the selection of the LAT events belonging 
to the TRANSIENT class, re-constructed under current 
the Instrument Response Functions (IRF) : Pass6_V3. 
Following FSSC suggestions, the selected energy range 
spanned from 100 MeV to 100 GeV. The data below 100 
MeV causes large uncertainties due to rapidly changing 
effective area with energy as well as ambiguity in the 
instrument response. We also set the maximum zenith 
angle value of 105° as suggested in the Fermi data pro- 
cessing "Cicerone 7 ". Finally, the time range for GRB 
data was set to be - 50 s after the burst trigger time 
(TO), almost the same as described in Akerlof et al. (A10, 
All). 



http: / /fermi. gsfc.nasa.gov/ssc/data/analysis/ 



7 http: / /fermi. gsfc.nasa.gov / ssc / data/analysis / documentation/C : 
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Fig. 2.— GBM and LAT light curves for GRB 081006A, in the 
order of increasing energy from top to bottom. The upper plot is 
the background subtracted sum of counts in the two Nal detectors, 
(NO + N3), in the 8 keV to 800 keV energy range with a 0.5 s bin 
size. The middle plot is the background subtracted counts in the 
BGO (B0) detector in the 280 keV to 4 MeV energy range with 
a 0.5 s bin size. The lower plot shows the number of LAT events 
which passed the TRANSIENT event selection above 100 MeV, 
bin step is 1 s. The open diamonds represent the energy for each 
photon. 
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The localization of GRB 081006A was also obtained 
from the TS map. The best fit position from the 
gtfindsrc procedure is RA = 135.53° and DEC = -61.65° 
(J2000). Figure 3 shows the error contours around the 
fitted position with 68%, 90% and 99% statistical error 
radii of 0.4°, 0.6° and 0.9°, respectively. This position 
is 0.55° from the location derived from the matched fil- 
ter technique. Thus the locations derived from the two 
methods are consistent to within 1 a error. 
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Fig. 3.— LAT localization for GRB 081006A with best fit po- 
sition RA = 135.53° and DEC = -61.65° (J2000). The contours 
around the estimated position shows a 68%, 90% and 99% statis- 
tical error radius of 0.4°, 0.6° and 0.9°, respectively. 



The likelihood analysis indicates a detection with a 
TS value of 45 assuming an isotropic background model. 
Considering that GRB 081006A is close to the Galactic 
plane (b ~ —0.9°), we must be careful that the high 
energy photon background does not have a local hot 
spot. However, a counts map with data taken three 
weeks prior to the burst demonstrated that the GRB 
position is safely >10° away from the bright region near 
the Galactic plane. We also estimated the count rate 
using the pre-burst data and re-scaled to the 50 s time 
range with TRANSIENT class events. Both these meth- 
ods show that the background is not more than 3 photons 
in the - 50 s time range. The TS detection value of 45 
for GRB 081006A is close to the significance of GRB 
080825C (Abdo et al. 2009d), GRB 081024B (Abdo ct 
al. 2010b) and GRB 090217A (Ackermann et al. 2010). 

In order to estimate the probability of such an occur- 
rence by chance alone, we performed identical searches 
using the likelihood analysis method on 16088 random 
LAT fields that were selected similarly to the procedure 
described in A10. Considering that the likelihood analy- 
sis is extremely time consuming, about 15 minutes on a 
PC, we only applied the likelihood analysis to the 2192 
random fields (out of the 16088) that contained more 
than three photons of class 2 or 3 within 16° radius re- 
gion. If there are less than three class 2 or 3 photons are 
within the region, it is almost impossible to detect a re- 
liable GRB signal. Figure 4 (upper panel) shows the cu- 
mulative distributions of the TS value for these 2192 ran- 



dom fields and the 27 (out of the 74) GRB fields that pass 
the same criteria. A K-S test indicates that the two dis- 
tributions are consistent if we exclude the outlier, GRB 
081006A. The random fields analysis shows that none 
have a TS value exceeding the candidate event, GRB 
081006A. The bottom panel of Figure 4 shows the cumu- 
lative distributions of the matched weight value (C^2 w i) 
derived from the statistical technique described in A10. 

In view of the heavy computational demands of the 
maximum likelihood method and the meager number 
of GRBs associated with high energy LAT photons, we 
decided to make some specific comparisons with the 
matched weight technique which is computationally at 
least a thousand times faster to compute. The distri- 
butions shown in Figure 4 indicate that the number of 
random fields that need to be examined to find a single 
event at the significance level of GRB 081006A is compu- 
tationally daunting. Without embarking on such a brutal 
course the most conservative comparison of the statisti- 
cal power of the two methods is obtained by noting that 
random probablity of obtaining a TS value greater than 
45 is less than 0.0123 = 27/2192 and the random proba- 
bility of obtaining a matched weight value greater than 
652 is less than 0.0043 = 74/17200. Thus, the matched 
filter technique is at least three times more effective in 
rejecting background. A more realistic estimate can be 
obtained by assuming exponential decreases of the cumu- 
lative distributions with the relevant parameters. These 
estimated trends are shown by the cyan lines in the two 
graphs. This would suggest that the matched filter re- 
jects background for similar events at a level at least 
50 times better than the maximum likelihood approach. 
Thus the false positive probability for this GRB identifi- 
cation is of the order of 10~ 5 . 

With the localization obtained from the matched filter 
technique, the likelihood analysis confirmed the identi- 
fication of faint LAT GRB 081006A. In several ways, 
these two methods are complimentary. The matched fil- 
ter is computationally efficient for identifying weak sig- 
nals with predetermined characteristics while the maxi- 
mum likelihood method works best when the combinato- 
rial background is relatively easy to estimate. The overall 
success of this approach is illustrated by GRB 081006A 
as well as the discoveries reported in previous papers. 

3.3. Joint spectral fitting with GBM and LAT data 

Joint spectral fitting with the GBM and LAT data was 
performed for the time range, -1 s to 8 s, after the trigger. 
Due to the faintness of the burst, time-resolved spectral 
fitting is not appropriate. RMFIT (version 3.3) software 
was used for spectral fitting with binned GBM time- 
tagged event (TTE) data. The LAT data is re-binned 
to to match the format required by RMFIT. 

Figure 5 shows the simultaneous fitting of the time- 
averaged count spectra from the GBM and LAT data. 
The Band function model (Band et al. 1993) fit to the 
GBM and LAT data gives a reasonable good fit with 
amplitude A = 2.1±g;| 10 ~ 3 (ph cm" 2 s" 1 kcV" 1 ), E p 



817.0+^40 (kcV), a = 0.78; 



0.24 



and /3 



n O o +0.10 
^8 _o,14 



(X /dof = 464/379). The corresponding energy flux is 
2.18±0.22 10" 7 erg cm" 2 s" 1 in the 10 keV to 1 MeV 
range and 6.06±0.61 10~ 7 erg cm~ 2 s _1 in the 10 keV 
to 1 GcV range. The E p of GRB 081006A is not well 
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constrained, as already noticed by van der Horst (2008) 
but the Band function model fitting is clearly better than 
a single power-law model or a power-law with exponen- 
tial cutoff model. Given that some GRBs have more 
complex spectral behavior than described by the Band 
function (e.g. Zhang et al. 2011), we also tried a two 
component Band function plus power-law, and also the 
Band function plus thermal component but the fit was 
not improved. 
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Fig. 4. — Upper panel : Complements of the cumulative distri- 
butions for TS values of 27 GBM GRBs (blue) and 2,192 similar 
fields obtained at random times (green). Bottom panel : The cu- 
mulative distributions for the matched weight values (CE^i) of 74 
GBM GRBs (blue), 1,731 random fields obtained nearly simultane- 
ously with the GBM data (red), and 17,200 random fields obtained 
at random times (green) . The dashed lines (cyan) in the two panels 
are extrapolations based on fits described in the text. (The dotted 
lines indicate estimated statistical uncertainties.) 



Faint high-energy gamma-ray photon emission of GRB 081006A from Fermi observations 
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Fig. 5. — Spectral distribution of time-integrated (-1 s to 8 s) 
GRB 081006A photons obtained from GBM and LAT data. The 
spectrum is well fit by a Band function model spanning ~5 decades 
of energy. The data from each instrument is indicated by color: red 
= NAI 0, blue = NAI 3, pink = BGO and green = LAT. The 
dashed lines depict the fit obtained with the Band function model 
convolved with the four different instrument response functions. 
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4. DISCUSSION 

GRB 081006A is detected with a fluence similar to 
three other GRBs, namely GRB 080825C, 081024B and 
090217A. Among these events, GRB 081024B is a short 
duration burst while GRB 080825C and GRB 090217A 
are long duration. We thus compare the properties of 
GRB 08 1006 A, also a long duration burst, with GRB 
080825C and GRB 090217A. 

The overall spectral features of GRB 081006A arc quite 
similar to GRB 080825C and 090217A: they can be fit- 
ted by the Band function model over five decades energy 
range from the 8 KeV to 800 MeV. No cut-off feature in 
the LAT energy range is detected. This suggests that a 
single emission mechanism is responsible for the broad- 
band emission of both GBM and LAT, e.g. synchrotron 
or jitter emission from the internal shock (Medvedev 
2000). 

The delayed onset of high energy LAT emission (>100 
MeV), observed for many other LAT-detected GRBs ( 
e.g. GRB 080825C, Abdo et al. 2009d; GRB 080916C, 
Abdo et al. 2009b; GRB 090510, Abdo et al. 2009c; 
GRB 090902B Abdo et al. 2009a), also shows a hint 
in GRB 081006A. The initial peak of LAT high energy 
photons with GRB 081006A, though correlated with the 
GBM peak, is delayed for about 2 s after the low energy 
GBM trigger. The delayed emission of high energy pho- 
tons from GRB 081006A is not as pronounced as found 
for GRB 080825C. In the later case, the onset of LAT 
emission is occurs at about 3 s and the first LAT peak 
is coincident with the second GBM peak (Abdo et al. 
2009d), but it is different from GRB 090217A for which 
there was no perceptible delay (Ackermann et al. 2010). 

Many GRBs ( e.g. GRB 080825C, Abdo et al. 2009d; 
GRB 080916C, Abdo et al. 2009b; GRB 090510, Abdo 
et al. 2009c; GRB 090902B Abdo et al. 2009a) also have 
long-lived high energy photon emission detected by the 
LAT even hours after the burst. The high energy photon 
emission of GRB 081006A clearly lasts longer (~ 40 s in 
LAT data) than the low energy emission (GBM Tgo ~ 7 
s). The two highest energy photons, both with E <~ 780 
MeV, are detected at 13.25 s and 26.5 s after the trigger 
as shown in Table 2. This is also similar to GRB 080825C 
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for which the LAT emission lasted slightly longer (up to 
T +35 s) than for the GBM (T 90 = 27 s). GRB 090217A 
does not show such longer higher energy emission. 

5. SUMMARY 

We have demonstrated the association of GRB 
081006A with high energy photon emission by apply- 
ing the matched filter technique to the Fermi /LAT 
data. The false positive probability is definitely less 
than 4 x 10~ 3 and probably much smaller. This event 
is found to be correlated with at least 13 high energy 
photons detected by the LAT instrument. A maximum 
likelihood analysis reveals a similar confidence level with 
a TS value of 45. Comparing the temporal and spec- 
tral properties with the other two long duration GRBs 
with similar fluences, GRB 080825C and 090217A, wc 
find GRB 081006A is closer to GRB 080825C. The delay 
and long emission duration for high energy photons are 
seen in GRB 081006A, similar to GRB 080825C, but not 
with GRB 0902 17A. These properties can be examined 
in more detail as the Fermi mission continues to obtain a 
larger sample of GRBs, especially faint ones such as GRB 
081006A, 080825C and 090217A. As demonstrated here, 
the matched filter technique is considerably more sensi- 
tive than the maximum likelihood analysis to find these 
fainter events. In the particular case of GRB 081006A, 
the background rejection is approximately 50 times bet- 
ter, making it a far better search tool for these rare 
events. As we have shown with several original iden- 
tifications of faint LAT GRBs, the two methods are best 
employed sequentially to first find the events and, sec- 
ondly, to determine the event characteristics. In this 
unique period when both Swift and Fermi are opera- 
tional, there is some urgency to explore the surprising 
dearth of faint LAT GRBs as fully as possible. 
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